Detection of high energy cosmic rays with the 
resonant gravitational wave detector 
NAUTILUS and EXPLORER 



P. Astone a , D. Babusci b , M. Bassan c ' d , P. Bonifazi a ' e , 
G. Cavallari h , E. Coccia c ' d , S. D' Antonio , V. Fafone c ' d , 

G. Giordano 13 , C. Ligi b , A. Marini b , G. Mazzitelli b , 
Y. Minenkov f , I. Modena f , G. Modestino b , A. Moleti c ' d , 
G. V. Pallottino a ' g , G. Pizzella b ' c , L. Quintieri b , A. Rocchi c , 
F. Ronga b '*, R. Terenzi c ' c , M. Visco c ' e 

a INFN (Istituto Nazionale di Fisica Nucleare) sezione di Roma 1 ,1 00185 Roma, 

Raly 

h INFN Laboratori Nazionali di Frascati, I 00044 Frascati, Raly 
c INFN Sezione di Roma II, I 00133 Roma, Italy 
d Dipartimento di Fisica, Universitd di Roma "Tor Vergata" , I 00133 Roma, Italy 
°Istituto di Fisica dello Spazio Interplanetario INAF, I 00133 Roma, Italy 
{ INFN Laboratori Nazionali del Gran Sasso, I 61010 Assergi, Italy 
g Dipartimento di Fisica, Universitd di Roma "La Sapienza" , I 00185 Roma, Italy 
h CERN, CH 1211 Geneva, Switzerland 



Abstract 

The cryogenic resonant gravitational wave detectors NAUTILUS and EXPLORER, 
made of an aluminum alloy bar, can detect cosmic ray showers. At temperatures 
above 1 K, when the material is in the normal conducting state, the measured 
signals are in good agreement with the values expected based on the cosmic rays 
data and on the thermo-acoustic model. When NAUTILUS was operated at the 
temperature of 0.14 K, in superconductive state, large signals produced by cosmic 
ray interactions, more energetic than expected, were recorded. The NAUTILUS data 
in this case are in agreement with the measurements done by a dedicated experiment 
on a particle beam. The biggest recorded event was in EXPLORER and excited the 
first longitudinal mode to a vibrational energy of ~ 670 K, corresponding to ~ 360 
TeV absorbed in the bar. Cosmic rays can be an important background in future 
acoustic detectors of improved sensitivity. At present, they represent a useful tool 
to verify the gravitational wave antenna performance. 
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1 Introduction 

Cosmic ray showers can excite mechanical vibrations in a metallic cylinder 
at its resonance frequencies and can provide an accidental background for 
experiments searching gravitational waves (gw): this possibility was suggested 
many years ago and a first search, ending with a null result, was carried out 
with room temperature Weber type resonant bar detectors [1]. 

More recently, the cryogenic resonant gw detector NAUTILUS has been equipped 
with shower detectors and the interaction of cosmic ray with the antenna has 
been studied in detail. 

The first detection of cosmic ray signals in a gw detector took place in 1998 
in NAUTILUS. During this run many events of very large amplitude were de- 
tected. This unexpected result suggested in 2002 the construction of a cosmic 
ray detector even for the EXPLORER detector. 

In section 2 we briefly recall the main features of the Thermo- Acoustic Model 
(TAM), that successfully describes the interaction between a solid elastic res- 
onator and a charged particle, or a beam of such particles; some of these 
features are extended to the regime of superconducting metal for the elastic 
resonator. 

In section 3 we describe the NAUTILUS and EXPLORER cosmic ray detec- 
tors and we specialize the TAM model to the interaction with cosmic rays 
computing the expected event rates. 

In section 4 we describe the results of coincidence measurements between the 
output of each gw antenna and its respective cosmic ray monitor, in different 
periods of data taking: for NAUTILUS during the year 1998, with the antenna 
in superconductive state, and then in the years from 2003 to 2006, while for 
EXPLORER in the period from 2003 to 2006. The data are interpreted with 
the help of some results obtained by the RAP experiment. Finally, some con- 
clusions of this extended analysis are drawn: we show the good agreement 
of our data with the TAM predictions and the consistency of data taken by 
two detectors with various different experimental setups (temperature, band- 
width, readout and acquisition hardware and software). As a central result, 
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the relevance of the conducting state of the antenna material on the strength 
of the interaction is proven. 

We have then learned that cosmic rays can also provide a good calibration 
source for resonant gw detectors, as they very closely mimic the signal ex- 
pected by short bursts of gravitational waves, i.e. the tidal excitation of the 
longitudinal modes. 



2 The thermo-acoustic model and its experimental validation with 
particle beams 

The interaction of energetic charged particles with a normal mode of an ex- 
tended elastic cylinder has been extensively studied over the years, both on 
the theoretical and on the experimental aspect. 

The first experiments aiming to detect mechanical oscillations in metallic tar- 
gets due to impinging elementary particles were carried out by Beron and 
Hofstander as early as in 1969 [2,3]. A few years later, Strini et al. [4] carried 
out an experiment with a small metallic cylinder and measured the cylin- 
der oscillations. The authors compared the data against the TAM (Thermo 
Acoustic Model) in which the longitudinal vibrations are originated from the 
local thermal expansion caused by the warming up due to the energy lost 
by the particles crossing the material. In particular, the vibration amplitude 
is directly proportional to the ratio of two thermophysical parameters of the 
material, namely the thermal expansion coefficient and the specific heat at 
constant volume. The ratio of these two quantities appears in the definition 
of the Gruneisen parameter 7. It turns out that while the two thermophysical 
parameters vary with temperature, 7 practically does not, provided the tem- 
perature is above the material superconducting (s) state critical temperature. 

Detailed calculations, successively refined by several authors [5,6,9,7,8] agree 
in predicting, for the excitation energy E of the fundamental vibrational mode 
of an aluminum cylindrical bar, the following equation: 

4 7 2 dW 2 nz sin[(7rl cos(9 )/2L} 2 
~ 97rpLvi [ dx ) 11 L ' nRcos(0 o )/L 1 [ ' 

where L is the bar length, R the bar radius, l the length of the particle 
track inside the bar, z Q the distance of the track mid point from one end of 
the bar, O the angle between the particle track and the axis of the bar, ^ 
the energy loss of the particle in the bar, p the density, v the longitudinal 
sound velocity in the material . This relation is valid for the material normal- 
conducting (n) state and some authors (see ref. [5,6]) have extended the model 
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to a super-conducting (s) resonator, according to a scenario in which the vi- 
bration amplitude is due to two pressure sources, one due to s — n transitions 
in small regions centered around the interacting particle tracks and the other 
due to thermal effects in these regions now in the n state. It is important to 
note, at this point, that a gw bar antenna, used as particle detector, has char- 
acteristics very different from the usual particle detectors which are sensitive 
only to ionization losses [9] [10]: indeed an acoustic resonator can be seen as a 
zero threshold calorimeter, sensitive to a vast range of energy loss processes. 



As anticipated in the introduction, the first detection of signals in a gw de- 
tector output due to cosmic ray events, took place in 1998. The NAUTILUS 
detector, a bar made of the aluminum alloy Al 5056 was operated at a thermo- 
dynamic temperature T = 0.14 K [11], i.e. below the s transition temperature 
T c ~ 0.9K. During this run, many events of unexpectedly large amplitude 
were detected. This result suggested an anomaly either in the model or in 
the cosmic ray interactions[12]. However the observation was not confirmed in 
the 2001 run with NAUTILUS at T = 1.5 K [13] and therefore we made the 
hypothesis that the unexpected behavior be due to the superconducting state 
of the material. An experiment (RAP) [14] was then planned at the INFN 
Frascati National Laboratory to study the vibration amplitude of a small Al 
5056 bar caused by the hits of a 510 MeV electron beam. The experiment 
was also motivated by the lack of complete knowledge of the thermophysi- 
cal parameters of the alloy Al 5056 at low and ultra-low temperatures. We 
summarize here the main results obtained by the experiment: 

(a ) in Al 5056 at T ~ AK [14] RAP measured a ratio, a n = 1.15, between 
the measured vibration amplitudes and the expectations based on the thermo- 
acoustic model of eq.l using the TAM parameters known for pure aluminum; 

(b) the experimental verification, made on a pure niobium bar, that the am- 
plitude depends on the material conduction state [15] and 

(c) for Al 5056 at T ~ 0.5 K a ratio, a s = 3.7, between the measured ampli- 
tudes at T ~ 0.5K and T ~ 1.5K [16] has been measured. 

While a n can be considered a small correction due to our inexact knowledge of 
relevant thermophysical parameters of our material, a value of a s so different 
from unity indicates that more complex interactions, as mentioned above, take 
place in the superconducting alloy. 
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3 The cosmic rays detectors of NAUTILUS and EXPLORER, and 
their expected rates 



The gw detector NAUTILUS [17] is located in Frascati (Italy) National Labo- 
ratories of INFN, at about 200 meters above sea level. It is equipped with a 
cosmic ray detection telescope made of seven layers of gas detectors (streamer 
tubes) for a total of 116 counters [18]. Three superimposed layers, each with 
an area of 36 m 2 , are located above the cryostat. Four superimposed layers are 
below the cryostat, each with area of 16.5 m 2 . The signal from each counter 
is digitized to measure the charge that is proportional to the number of parti- 
cles. The detector is capable of measuring particle densities up to 1000 vart ^ es 
without large saturation effects. During normal runs only showers are detected, 
with a typical threshold of the order of 2 par ^ les in the lower detectors. 

Single particles, cosmic rays muons, are collected every day in short runs, with 
the aim of calibrating the detectors. The systematic error on the absolute 
number of particles crossing the apparatus is of the order of 25%. The cosmic 
ray and the antenna data acquisition systems are independent. A GPS clock 
is used to synchronize the antenna and cosmic ray acquisitions. The time 
resolution is limited to 0.2 ms by the antenna ADC 5 kHz sampling. 



The gw detector EXPLORER [20] is located in CERN (Geneva-CH) at about 
430 meters above sea level. Scintillators counters were installed at EXPLORER 
in 2002, using scrap equipment recovered after the LEP shutdown. Above the 
cryostat there is a single layer of 11 scintillators for a total area of 9.9 m 2 . 
Below the cryostat there are two layers of 4 counters each, with a total area 
of 6.3 m 2 . Each scintillator is seen by two photomultipliers (56AVP). The 
signals from the anode and from the last dynode of each photomultiplier are 
digitized to measure the total charge. No large saturation effects occur in the 
dynodes up to particle densities of the order of 2000 pctr ^ es - The typical trigger 
threshold during normal run is of the order of ^ P artt f es j n the lower detectors. 

The detector is calibrated using cosmic ray muons as in NAUTILUS. Correc- 
tions are applied to take into account the scintillators attenuation lengths and 
the conversion of photons inside the scintillators. The cosmic ray and the an- 
tenna data acquisition systems are independent, as in NAUTILUS. Similarly, 
a GPS clock is used to synchronize the acquisitions. 



In order to compare the shower particle densities measured by the scintillators 
and by the streamer tubes, we have installed two scintillators, equal to the 
ones used in the top layer of EXPLORER, above the NAUTILUS cryostat. 
The scintillators measured numbers of particles +20% larger than streamer 
tubes for showers particle densities around 400 P artl f es , 
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Fig. 1. The integral rate of events computed using Cocconi's relation (eq.2) (con- 
tinuous line), the rate measured by the EXPLORER scintillators (squares) and by 
the NAUTILUS streamer tubes (circles). The open symbols refer to the detectors 
below the cryostat, the filled symbols are for the detectors above. The EXPLORER 
data are corrected for the 20% systematic difference in the number of particle mea- 
sured with scintillators with respect to streamer tubes. The difference in the rates 
of the detectors above the cryostats is due to the different altitudes of Frascati and 
CERN and to the presence of a roof in the EXPLORER building. The rates in the 
bottom detectors are also influenced by the different distribution of iron, copper 
and aluminum in the cryostats. 

Most of the high energy events are due to electromagnetic showers. The rate 
of electromagnetic air showers (EAS) is computed starting from the empirical 
relation due to G. Cocconi [21] 

H(> A) = A;A~ A events/day (2) 



where A is the density of secondaries in an EAS, measured in units of number of 
charged particles per square meter, A = 1.32 + 0.038 ln(A) and k = 3.54 x 10 4 . 
This relation holds at sea level and in absence of absorbing material. The 
NAUTILUS antenna is located inside a building with a very small amount of 
matter in the roof, while EXPLORER is in a normal building with concrete 
roof. Concrete has 50 MeV critical energy, to be compared to a critical energy 
88 MeV for air; therefore we expect in the EXPLORER detector an increase 
of the electromagnetic showers particle density due to the different critical 
energies and to its above sea level higher location. 

The integral distribution of the showers measured with both EXPLORER 
and NAUTILUS cosmic ray detectors is shown in Fig.l. In this figure we 
have (arbitrarily) taken NAUTILUS as reference: therefore we have shifted 
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the EXPLORER points in order to correct the systematic 20% error between 
streamer and scintillators. 

The NAUTILUS data of the detector above the cryostat are in agreement with 
the prediction of Eq.2 within the 25% systematic error given by the particle 
density measurement. The difference between NAUTILUS and EXPLORER 
and between detectors over and under the cryostats are due to several, already 
mentioned effects: the differences in altitude of the experimental locations 
(230 m), the presence of a concrete roof in the EXPLORER building and the 
materials in the cryostats. 

The signal expected in a gw detector like NAUTILUS, a bar 3 m in length 
and 0.6 m diameter, as a consequence of the interaction of a particle releasing 
an energy (W in GeV units) is [11,12,13] , according to relation (1): 

E ~ Z|l • 10~ 9 W 2 a 2 [K[ (3) 

where the bar oscillation energy E is expressed, as usual in the antenna jargon, 
in kelvin units (IK = 1.38 • 10~ 23 J), the numerical constant is the value 
computed using the linear expansion coefficient and the specific heat of pure 
aluminum at 4 K and a takes, as described in the previous section, either value 
a n = 1.15 above the s transition temperature or a s = 3.7 for superconductive 
Al 5056. The constant 7.64 • 10 -9 applies if the energy is released in the bar 
center. If the energy is uniformly distributed along the bar, as in the case of 
EAS showers, this value is reduced by a factor 2. 

The cosmic ray event rate in NAUTILUS has been evaluated considering three 
different event categories: pure electromagnetic showers, showers produced by 
muons and showers produced by hadrons in the bar. We use Eq.3 with the 
correction a n = 1.15 for the response of an aluminum Al 5056 bar in the 
normal state. 

The rate of the EAS and the energy deposited by an EAS has been computed 
starting from Eqs. 2,3 with the following assumptions: 

1) No particle absorbed (all particles go through the bar): indeed the radia- 
tion length in the bar is small compared to the total radiation length in the 
atmosphere. 

2) The energy loss for a single particle is computed assuming ionization energy 
losses for electrons having the aluminum critical energy. 

3) We used the showers angular distribution as reported in[22]. 

4) We neglected the contribution of hadrons that could be present in the core 
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of the showers. 



Under the previous assumptions and using the density A of secondaries we 
obtain [11,12,13]: 

E = A 2 4.7 1(T 1( V [K[ (4) 



The production of the showers due to muon and hadrons was computed using 
the GEANT package [25], developed at CERN, to simulate NAUTILUS and the 
CORSIKA[23] Montecarlo, as input to GEANT, to simulate the effect of the 
hadrons produced by the cosmic ray interactions in the atmosphere, assuming 
a cosmic ray "light" composition. The Montecarlo simulation reflects 1 year 
of data taking. 

The results are shown in Table 1. The rate of the the events scales as W~ - 9 . 
This is because the cosmic ray integral spectrum is well described by a power 
law with (5 ~ 1.7 for cosmic ray primaries up to the so called "knee" 

at W pr i mar y = 10 15 eU and (3 ~ 2 at higher energies. The energy in the first 
longitudinal mode E (first column of Table 1) is proportional to the square of 
the absorbed energy W. 
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Table 1 

Estimated rate (events/day) of antenna excitations due to cosmic rays in NAU- 
TILUS as a function of the vibrational energy of the longitudinal fundamental mode 
that such events can produce. The value at E = 0.1AT is obtained extrapolating 
from the lowest energy values. The values in the second column are the energies ab- 
sorbed by the bar computed from Eq. 4, with the assumption of energy uniformly 
distributed, and a n = 1.15. 

There is quite a large uncertainty in the estimation of the high energy event 
rate. This is due to uncertainties both in the cosmic ray composition and in 
the models of hadronic interactions at high energies. We have performed a 
check of the hadron flux at sea level used in our simulation with the direct 
hadron flux measured by the EAS-TOP experiment [24], properly scaling their 
results for the different altitude. We have found that at 1 TeV the EAS- 
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Fig. 2. NAUTILUS 1998. Distribution of the measured E exp , with A > 50 par ^ les . 
The upper graph shows all data. In the lower graph a cut for T e jj < 5 mK was 
applied. The negative values are due to the subtraction of the noise energy (T e fj in 
the antenna jargon). 



TOP measurement gives a flux roughly +35% higher that the one used in 
our simulation. This gives an idea of the uncertainty in the simulation of the 
hadronic effect. 



Another uncertainty is due to the EAS rate, modified by the presence of ma- 
terials, as shown in Fig.l. Comparing the rates of the EXPLORER and NAU- 
TILUS lower detectors and extrapolating Fig.l at higher densities we have 
estimated that EXPLORER should have an excess of events respect to NAU- 
TILUS of a factor 2.8 for energies larger than 0.1 K. We underline that, due 
to the large uncertainties involved, the expected absolute rate of events pro- 
ducing signals in a gravitational wave bar has also a large uncertainty, of the 
order of that shown in Fig.l. These, however affect in the same manner both 
our antennas, so that the uncertainty on the relative rates of EXPLORER and 
NAUTILUS is much smaller, being only due to systematic errors in the cali- 
bration of the EAS detectors (~ 25%) and of the gravitational wave detectors 
(~ 10%). In the following, we shall use only the particle densities measured 
by the lower detectors as they are closer to the bar. 
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4 Antenna Signals Generated by Cosmic Rays 



4.1 NAUTILUS in 1998 



The ultra-cryogenic resonant-mass gravitational wave (gw) detector NAU- 
TILUS [19] operating since 1996 at the INFN Frascati Laboratory, consists 
of a 3 m 2300 kg Al 5056 alloy bar. The cryostat mainly consists of seven con- 
centric layers: three steel vessels, two thin aluminum plus three thick copper 
thermal shields. During the run of 1998 it was cooled at 140 mK. The quan- 
tity that is observed (the "gw antenna output") is the vibrational amplitude 
of its first longitudinal mode of oscillation. This is converted by means of an 
electromechanical resonant transducer into an electrical signal which is am- 
plified by a dc-SQUID. The bar and the resonant transducer form a coupled 
oscillator system, with two resonant modes, whose frequencies were, in 1998 
/_ = 906.40 Hz and /+ = 921.95 Hz. 

The data regarding the vibrational energy of the NAUTILUS gw antenna 
were recorded with a sampling time of 4.54 ms and processed with the delta- 
matched filter [26] optimized to detect impulsive signals. In a previous paper 
[11] we reported the results of a search for correlations between the NAUTILUS 
data and the data of the EAS detector, when for the first time acoustic signals 
generated by EAS were measured. In a further investigation [12] , we found very 
large NAUTILUS signals at a rate much greater than expected. Now we know 
that, since the bar temperature was about 0.14 K, the value a s = 3.7 must be 
used in Eq. 4 to compute the expected response. 

The correlation between the small signals detected by NAUTILUS and the 
impinging EAS has been described in detail in [12]. The main points of this 
procedure are: 

(1) We consider stretches of the filtered NAUTILUS antenna data corre- 
sponding to EAS with density A , with a lower threshold A > 50 par ^ 2 des . 

(2) For each stretch we calculate the average energy E, in a time interval 
±227 ms around the EAS arrival time, subtracting the value due to the 
noise energy (T e ff in the antenna jargon). The time interval is chosen to 
take into account the expected shape of the offline filtered signal. 

(3) With this averaging procedure we avoid the problem of taking either a 
maximum or a minimum value, which may be due to noise and, when due 
to signals, might not be exactly in phase among the various stretches. By 
doing so we get average values E. In order to convert the value E into the 
energy at the maximum E exp , we multiply E by a factor 4.1, as found, 
with a statistical dispersion of a few percent, by numerically averaging 
the data sample of big events where the signal is much larger than the 
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Fig. 3. Averages of signals with energy E exp < 0.1 K, grouping data in ranges of 
particle density A. Filled circles NAUTILUS at T = 0.14 K, open circles NAU- 
TILUS at T = 3 K, filled squares EXPLORER at T = 3 K. The data gathered 
at T = 0.14 K are almost one order of magnitude larger than those collected at 

r = 3 k. 

background, so that noise effects can be neglected. 
(4) We obtain 26425 stretches of filtered data in coincidence with EAS. The 
NAUTILUS average noise level during this run was T e ff ~ 10 mK, 
while in the NAUTILUS 2003-2006 run the noise was T eff ~ 4 mK. In 
order to perform a more meaningful comparison of these data with the 
NAUTILUS 2003-2006 run we have considered only those stretches with 
T e ff < 5 mK. In this way the number of useful stretches of filtered data 
reduces to 8904, in the period October 1998 to December 1998, for a total 
live-time of 27.3 days. 

In order to verify the TAM model, we eliminate large signals with energy 
Eexp > 100 mK and we bin the remaining in five ranges according to the 
particle density A, measured by the streamer tubes under the cryostat with 
an upper cut to A = \QQQ P artl ^ les ^ Q avoid the saturation effects in the cosmic 
ray detectors. 

The plot of excitation energy E exp vs particle density A is shown in Fig. 3. In 
this figure we show both the measurements with NAUTILUS at 140mK and 
2.6K, as well as EXPLORER at a temperature of about 3 K (see discussion 
in the following sections). We clearly see a difference of almost an order of 
magnitude between the measurements taken with aluminum in the (s) state 
and those in (n) conduction state. 
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Fig. 4. NAUTILUS 1998: small signal analysis. Signals with energy E exp <0.1K. 
Measured average energy < E exp > of signals with their standard deviations in 
K units versus signal energy as calculated with Eq. 4, using the particle densities 
measured by the streamer tubes below the cryostat and with a s = 3.7. The straight 
line is a least square fit through the origin, the slope is pO = 0.82 ± 0.17. The fit is 
in good agreement with the predictions, given the systematic errors on the number 
of particles (~ 50% ) and on the antenna calibration (~ 10%) 

We apply Eq.4 to evaluate the expected response according to the TAM model 
with using a s = 3.7 . The results are shown in Fig. 4: the agreement between 
the measurements and the values calculated with the TAM is reasonably good. 
The slope of the linear fit passing through the origin is 0.82 ± 0.17 compatible 
with one, and this confirms that the model works well also in the (s) state. The 
systematic uncertainty in the number of particles is ~ 25%, therefore, as E 
scales quadratically with A (see Eq. 4) we expect an error in the determination 
of the slope of ~ 50%, plus a smaller (~ 10%) uncertainty due to the antenna 
calibration . 

We have also measured the rate of events producing signals in the bar. To this 
purpose we have used the energy distribution of Fig.2 with noise T e ff < 5 mK. 
This energy distribution is the convolution of the cosmic ray signals with the 
noise due to random fluctuations of the background. We have computed the 
background energy distribution by shifting the time of the cosmic rays (20 
different values with 2 s intervals for each cosmic ray). The event rate per 
day after the unfolding of the background distribution is shown in Fig. 5. The 
agreement with the predictions, computed from the figures of Table 1 modified 
by using the correct value of a s , appears very good (taking into account the 
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Fig. 5. NAUTILUS 1998. The integral distribution of the event rate after the 
background unfolding, compared with the expected distribution (continuous line). 
The prediction is computed using the data of Table 1 and using the appropriate 
value a s = 3.7. 

very large uncertainties in the expected rates). 
4.2 NA UTIL US in 2003-2006 

In 2003 some components in the readout of NAUTILUS antenna were up- 
graded, similarly to what has been done in EXPLORER the previous year 
[27]. This resulted in an increased bandwidth from the 1998 value of 0.4 Hz 
to about 9.6 Hz. The data were since recorded with a sampling time of 0.2 ms 
and processed with the delta-matched filter [26] . NAUTILUS is operated at 
about 3 K by cooling with a superffuid 4 He bath, so that a n = 1.15 is used 
in Eq.4 to estimate the signal expected according to the TAM. 

We have carried out the small signals analysis according the following proce- 
dure, similar to that outlined in of the previous section: 

(1) The data for the four years 2003-2006 were grouped together with a lower 
threshold A > 50 2^"^. 

(2) For each stretch of filtered data the average energy E is calculated using 
160 contiguous samples, corresponding to ±16 ms around the EAS arrival 
time. 

(3) The adjusting factor from the averages E to the maximum E exp is 2.3, as 
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Fig. 6. NAUTILUS in 2003-2006. Distribution of the measured antenna output 
Eexpi with A > 50 fo-rudes ^ rj-^g U pp e r graph reports all data. The lower graph is 
for T e ff < 5 mK. 

found by numerically averaging large events where noise contribution is 
negligible. The difference with respect to the previous value 4.1 is mainly 
due to the larger bandwidth of the detector and to the different time 
window. 

(4) We obtain in total 88263 stretches of filtered data, corresponding to EAS, 
and 82330 with T e ff < 5 mK. The live-time of the data after the cut 
T e ff < 5 mK is 1086 days. The distribution of E exp is shown in Fig.6 

With the same procedure used for the 1998 data, we consider 5 ranges ac- 
cording the multiplicity A > 50 vart ^ es and apply Eq.4 to compute the ex- 
pected values according to the TAM. The results are shown in Fig. 7. The 
ranges in this figure are smaller than the ones in Fig.4, as a consequence of 
the different values of a in superconductive and in normal state. The slope 
of a linear fit through the origin is 0.66 ± 0.09. This value should be com- 
pared with 0.64 ± 0.12, that we obtain by dividing for o? n the published value 
0.85±0.16[12]. 

The analysis of ref. [12] only concerned the data of NAUTILUS 2001: at that 
time NAUTILUS was still working in the narrow band mode with the old 
readout system. The excellent agreement demonstrates the stability of the 
apparatus and the accuracy of calibrations, even in largely different set-ups. 
Moreover we stress the good agreement between measured and calculated 
quantities, within the TAM. 
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Fig. 7. NAUTILUS 2003-2006: small signal analysis. Signals with energy 
Eexp < 0.1 if with their standard deviations compared to the signals calculated 
with Eq. 4 and a n using the particle densities measured by the streamer tubes un- 
der the cryostat. The straight line is a least square fit through the origin. The fit 
is in good agreement with the prediction taking into account the systematic errors 
from the measurement of the number of particles (~ 50% ) and from the antenna 
calibration ( ~ 10%) 

Figure 8 shows the integral distribution of the event rate after the background 
unfolding. The agreement with the expected distribution computed from Ta- 
ble 1 is again quite good, considering the large uncertainties in the expected 
rates. The largest event detected in NAUTILUS in the years 2003-2006 has 
an energy of E exp ~ 4.1 K in the first longitudinal mode corresponding to 
W ~ 28 TeV in the bar. 

13 EXPLORER m 2003-2006 

The EXPLORER detector has been in almost continuous operation at CERN 
since 1991, and it has undergone over the years several upgrades that progres- 
sively improved both its sensitivity and its operation duty cycle. EXPLORER 
has a bar similar to NAUTILUS, while the cryostat is sligthly different (three 
steel containers, one aluminum shield and a thin copper vessel). Data acqui- 
sition, readout and operation are very similar to NAUTILUS. The operating 
temperature is T ~ 2.6 K. A detailed description of the apparatus and its 
main features (including data taking and analysis) can be found in ref [20]. 
In 2001 EXPLORER has been upgraded with a new read-out allowing for the 
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Fig. 8. NAUTILUS 2003-2006 : The integral distribution of the event rate after the 
background unfolding, as in fig. 5, for the four years 2003-2006, compared with the 
expected distribution (continuous line). The prediction is computed using the data 
of Table 1. 
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Fig. 9. EXPLORER in 2003 2004 2005 2006. Distribution of the measured antenna 
output E exp , with A > 50 P ar ^ les , The upper graph for all data. The lower graph 
for T e jf < 5 mK. Note the very large event with energy bigger than 600 K in the 
upper graph. 
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Fig. 10. EXPLORER 2003-2006: small signal analysis. Signal energy of events with 
Eexp < 0.1 if with their standard deviations vs the energy calculated with Eq. 4 
and a n = 1.15, using the particle densities measured by the scintillators under the 
cryostat. The straight line is a least square fit through the origin. The fit is in good 
agreement with the prediction taking into account the systematic errors from the 
measurement of the number of particles (~ 50%) and from the antenna calibration 
(~ 10%). There is also a good agreement with NAUTILUS data. 

first time "wide band" operation of a gw bar detector [27]. Explorer operation 
was suspended in August 2002 due to a cryogenic failure. We took advantage 
of this stop to recondition the transducer and complete installation of the 
cosmic ray shower detector described in section 3. 

For the small signal analysis we have applied the same procedure used for 
NAUTILUS, with the following differences: 

• The adjusting factor from the averages E to the maximum E exp is 2.5, again 
determined by numerically averaging large events. 

• We obtain in total 431256 stretches of filtered data, corresponding to EAS, 
and 407064 with T e ff < 5 inK. The live-time of the data after the cut 
T e ff < 5 mK is 1022 days. The distribution of E exp is shown in Fig.9. 

Fig. 10 shows the small signal response analysis: as we did for NAUTILUS, 
we plot the predicted energy vs. the measured energies (both expressed in 
mK). The linear fit for the EXPLORER data is compatible with that in Fig. 7 
for NAUTILUS. This shows that we have in good control both the calibra- 
tion of the cosmic ray detectors and the calibration of the gravitational wave 
detectors. 



17 



EXPLORER 2003 2004 2005 2006 T=3 K 




Fig. 11. EXPLORER 2003-2006 : The integral distribution of the event rate after 
the background unfolding, as in Fig. 5, compared with the expected distribution 
(continuous line). The prediction is computed using Table 1 multiplied by a factor 
2.8 (see text). 



We also repeated the large event analysis and measured the rate of the large 
events. The results are in Figure 11, similar to Fig. 2 and 8. The event rate in 
EXPLORER is higher that in NAUTILUS. We expect a higher rate due the 
different altitude of Frascati and CERN and to the effect of the roof in the 
CERN building. The continuous line in Fig. 11 shows the predictions computed 
from Table 1 scaled by a factor 2.8 that accounts for the difference in the EAS 
rates as measured by the cosmic ray detectors and discussed in section 3 . 



The agreement between measurement and expectations is again quite good, 
considering the large uncertainties in the calculation of the predicted rates. It 
is important to note that acoustic gw detectors have no large signal limitations 
due to saturation effects and can detect very high energy events. 



Indeed the largest event detected up to now has an energy in the first longi- 
tudinal mode of ~ 670 K corresponding to ~ 360 TeV in the bar. The event 
occurred in EXPLORER on Nov 10 2006 9:40 UT. 
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Fig. 12. EXPLORER 2003-2006 : Time difference (seconds) between cosmic rays 
with A > 100 vari ^ es and the maximum of the filtered antenna signal, with a cut 
E > 36 T e ff. The fit with a gaussian , with parameters p0=peak, pl=mean, p2=<r 
and a constant background p3, gives a = 3.7ms. The value of the mean (-1±0.35 
ms) should be compared to the expected value of -0.6 ms due to the delay of the 
antenna electronic chain. 

5 Conclusions 



We have discussed the NAUTILUS and EXPLORER response to small signals 
(energy E < a few mK) and to large signals (from energy E > 20 mK up to 
events of E ~ 600i^). Table 2 shows the main results obtained. 



RUN 


T 


small signal 


expected 


large ev rate 


large ev rate 




(K) 


slope 


slope 


E > 0.1K 


predicted 


NA 1998 


0.14 


0.82 ±0.17 


1±0.5 


0.53 ±0.16 


0.24 


NA 2001 [13] 


1.5 


0.64 ±0.12 


1±0.5 


0.007 ±0.012 


0.03 


NA 2003-06 


~ 3 


0.66 ±0.09 


1±0.5 


0.034 ±0.007 


0.03 


EX 2003-06 


2.6 


0.6 ±0.07 


1±0.5 


0.078 ±0.01 


0.08 



Table 2 



Data summary. The column "small signal slope" shows the results of the fits in 
Fig. 4, 7 and 10. The predictions are based on RAP experiment [14] [15] [16]. The 
systematic error on the slope is ~ 50%. The errors reported in column three are 
statistical. The columns "large event rate" refer to the number of events per days 
having energy E > 0.1K. 
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For small signals we have found that the acoustic gw detector response well 
agrees with the predictions based on the Thermo-Acoustic Model, once the 
corrections to this model, provided by a dedicate experiment [14] [15] [16], are 
applied. 

The large signals are the kind of events that can represent a background noise 
for current gravitational wave detectors. We found the rate of large signals 
to be in good agreement with the predictions, given the large uncertainty in 
such predictions. We notice that in EXPLORER we have a higher rate than 
in NAUTILUS, and this is due to the the concrete roof and to the difference 
in altitude between Frascati and Geneva. 

We have shown that the unexpected large events detected in 1998 with NAU- 
TILUS at T=0.14 Kelvin were due to the superconductive state. 

Cosmic rays noise could become an important noise in higher sensitivity de- 
tectors , namely in superconductive state, and this noise should be taken 
into account in possible future detectors of improved sensitivity, both acoustic 
[29] [30] and interferometric. As shown in this paper, cosmic rays can provide 
an useful tool to have a continuous monitor and calibration of the acoustic 
gravitational wave detectors. As an example Fig. 12 shows the antenna time 
resolution measured using the cosmic ray showers collected during the EX- 
PLORER run. 

Very important is the fact that the observation of cosmic rays demonstrates 
that the detectors, both hardware and software, are indeed able to detect 
vibrations as small as 10~ 19 m. 
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